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Abstract 

We report on the results of detailed X-ray spectroscopy of the Fe-K region in the Seyfert 1 galaxy NGC 3783 
using the Chandra High Energy Grating Transmission Spectrometer [HETGS). There were five observations over 
an interval of ^ 125 days in 2001, each with an exposure time of ~ 170 ks. The combined data constitute 
the highest signal-to-noise Fe-K spectrum having the best velocity resolution in the Fe-K band to date (FWHM 
~ 1860 km s“^). The combined data show a resolved Fe Ka line core (FWHM = 1700jl39g with a center 

energy of 6.397±0.003 keV, consistent with an origin in neutral or lowly ionized Fe, located between the BLR and 
NLR, as found by Kaspi et al. (2002). We also find that excess flux around the base of the Fe Ka line core can be 
modeled with either a Compton scattering “shoulder” or an emission line (with about the same flux as the line 
core) from a relativistic accretion disk, having an inclination angle of 11° or less. This disk line model is as good 
as a Compton-shoulder model for the base of the Fe Kq line core. In the latter model, we measured the column 
density to be x 10^^ cm“^, which corresponds to a Thomson optical depth of ~ 0.60, so the line-emitting 

matter is not quite Compton-thick. An intrinsic width of ISOOlg^Q km s“^ FWHM is still required in this model. 
Moreover, more complicated scenarios involving both a Compton-shoulder and a disk line cannot be ruled out. 
We confirm an absorption feature due to He-like Fe (FWHM = 6405l|g7Q km s“^), found in previous studies. 
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1. INTRODUCTION 

At least part of the Fe Ka fluorescent emis¬ 
sion line in type I active galactic nuclei (AGNs) 
is believed to originate in a relativistic accretion 
disk around a black hole (e.g. see reviews by 
Fabian et al. 2000; Reynolds & Nowak 2003). 
The dominant peak energy of the Fe Ka line at 
~ 6.4 keV appears to be ubiquitous and this core 
of the line carries a substantial fraction of the to¬ 
tal line flux (e.g. Nandra et al. 1997; Sulentic 
et al. 1998; Lubihski & Zdziarski 2001; Weaver, 
Gelbord, & Yaqoob 2001; Yaqoob et al. 2002; 
Perola et al. 2002; Reeves 2003; Yaqoob & Pad- 
manabhan 2004). It has been traditional to asso¬ 
ciate such narrow Fe Ka lines with an origin in 
distant matter, at least several thousand gravi¬ 
tational radii from the putative black hole (e.g. 
the optical broad-line region (BLR), the puta¬ 
tive obscuring torus, or the optical narrow-line 
region (NLR)). However, Petrucci et al. (2002) 
recently reported a variable, narrow Fe Ka line in 


Mkn 841, supporting an accretion-disk origin (see 
also Longinotti et al. 2004). Moreover, rapidly 
variable, narrow Fe K line emission has been ob¬ 
served in the Seyfert I galaxy NGG 7314 (Yaqoob 
et al. 2003a). Thus, even narrow Fe Ka lines 
may have a significant contribution from the ac¬ 
cretion disk (Lee et al. 2002; Yaqoob et al. 2003a; 
Longinotti et al. 2004; Turner, Kramer, & Reeves 
2004). 

NGC 3783 is a fairly bright (F 2 _io keV ~ 
5 — 7.5 X 10“^^ ergs s“^), moderate luminosity 
(L2_io keV ~ 1 — 1.5 X 10"^^ ergs s“^), nearby 
{z = 0.00973) Seyfert 1 galaxy which has been 
well-studied in all wavebands. In the X-ray band, 
NGC 3783 was the target of the deepest observa¬ 
tion with the Chandra high energy grating trans¬ 
mission spectrometer (HETGS) for any Seyfert 
galaxy, during a campaign in 2001, obtaining a 
net exposure time of ^ 830 ks. The spectra from 
this campaign represent the highest spectral res¬ 
olution data with the best signal-to-noise for any 
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Table 1 

Core Fe K Line Chandra (HEG) Spectral Fitting Results 


Observation 



EW" 

FWHM'^ 



(keV) 


(eV) 

(km s“^) 



NGG 3783(1) 

6 491 +0 008 

D.UUi_Q QQg 

4 77-I-1.31 
+• ‘ ‘ -1.17 

7 c :+21 

'0-18 

9q9f»+1225 

z:oz,u_g95 

5.0 


(6.390 - 6.412) 

(3.22 - 6.58) 

(51 - 103) 

(990 - 4055) 

1.1 

NGG 3783(2) 

6 303+0 008 

U.0»0_g QQg 

5.42+)'39 

85l?? 

25451)?°° 

5.1 


(6.380 - 6.403) 

(3.82 - 7.33) 

(60- 115) 

(1365 - 4415) 

1.1 

NGG 3783(3) 

6 304+0 005 

O.OaU_g ggg 

5 IQ+i-^o 

8II)? 

1 205+°^° 
1ZUO_j205 

5.2 


(6.387 - 6.401) 

(3.77- 6.83) 

(59 - 107) 

(0 - 2540) 

1.1 

NGG 3783(4) 

6 305+0 005 

U.O»U_g QQg 

5.611);?® 

661)® 

1255l)??5 

7.4 


(6.389 - 6.402) 

(4.08 - 7.31) 

(48 - 86) 

(0 - 2305) 

1.6 

NGG 3783(5) 

6 401+° °°^ 

O.UUl_g QQg 

3.86l);« 

ri+15 

OJ^-13 

1 343+990 
iOUO_j34g 

6.1 


(6.392 - 6.411) 

(2.52 - 5.39) 

(33 - 71) 

(0-2715) 

1.3 

Total 

6 307+0 003 
u.o»( -0.003 

4 00+0-55 
u.au_g 52 

70l? 

17001"° 

5.8 


(6.393 - 6.401) 

(4.21 - 5.64) 

(60 - 80) 

(1180 - 2250) 

1.2 


Chandra HEG data, fitted with a power law plus Gaussian Fe Ka emission-line model in the 2-7 keV band (see 
§3). The purpose of these fits is primarily a comparison of the Fe Ka emission-line parameters (see Fig. 3, Fig. 4, 
and Fig. 6 for more complex fits). Absolute values of the line intensity are model-dependent but the peak energy 
and width of the Fe Ka line core are insensitive to the continuum and absorption model. The last row gives 
measurements from the summed spectrum (with ~ 830 ks exposure time), as derived in Yaqoob & Padmanabhan 
(2004). All fit parameters are quoted in the source rest frame. Line widths are intrinsic, since they have already 
been corrected for the instrumental line-response function. Statistical errors are for three interesting parameters 
and correspond to the 68% confidence level (AC-statistic = 3.51), whilst parentheses show the 90% confidence 
level ranges of the parameters (AC-statistic = 6.25). “ Gaussian line center energy. ^ Emission-line intensity 
in units of 10“® photons cm~^ s“^. Emission line equivalent width. Full width half maximum, rounded to 
5 km s“^. ® F is the estimated 2-10 keV observed flux in units of 10~^^ ergs cm“^ s“^ using the complex disk-line 
model in §6. The power-law continuum was extrapolated to 10 keV. L is the estimated 2-10 keV source-frame 
luminosity (using the 2-10 keV estimated flux), in units of 10^® ergs s“^, assuming Hq = 70 km s“^ Mpc“^, 
flmatter ~ 0.3 and Ha — 0.7. 


Seyfert galaxy available to date. The observing 
campaign was also designed to study variability 
and was broken up into five observations sepa¬ 
rated by various intervals ranging from days to 
months. The soft X-ray part of the HETGS data 
has been well-studied by several research groups 
and has produced a wealth of new information 
and insights into the photoionized outflow and 
its variability (Kaspi et al. 2001, 2002, Netzer 
et al. 2003, Krongold et al. 2003). These re¬ 
sults have been supplemented by studies using 
XMM-Newton (Behar et al. 2003; Reeves et al. 
2004), and by studies in the UV band (e.g. Gabel 
et al. 2003a,b) On the other hand, the high- 
spectral resolution HETGS data for the Fe Ka 


line in NGC 3783 from the extended campaign 
in 2001 have not yet been fully exploited. Kaspi 
et al. (2002) reported measurements of the nar¬ 
row core of the Fe Ka line and detection of a 
Gompton-scattering “shoulder” on the red-side of 
the core, indicating an origin for the line core in 
cold, optically-thick matter, far from the nucleus 
(beyond the optical BLR). Results of variability 
studies of the Fe Ka line from within the extended 
observation campaign in 2001 have not yet been 
reported. A detailed study of the Fe Ka line in 
NGG 3783 using XMM-Newton was presented by 
Reeves et al. (2004) who confirmed the Chandra 
measurements of the narrow core, and in addi¬ 
tion reported that a broad, relativistic disk line 
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was not required by the data (a similar conclu¬ 
sion was reached by Kaspi et al. (2002) for the 
Chandra HETGS data). 

In this paper we specifically study the Fe Ka 
line in detail in NGC 3783 using the Chandra 
HETGS data from the extended campaign in 
2001. We stndy the time-averaged spectrum, as 
well as the individual spectra from the five ob¬ 
servations from the campaign in order to address 
line variability. We find that a Compton shoul¬ 
der model is not a unique description of the com¬ 
plexity in the Fe Ka line. A relativistic disk 
line model provides as good a fit as the Compton 
shoulder model. In §2 we describe the data and 
observations. In §3 we discuss the results of spec¬ 
tral fitting to derive parameters for the Fe Ka line 
core, from the time-averaged spectrum as well as 
from the separate observations. In §4 we confirm 
the detection of an Fe He-like absorption feature. 
In §5 and §6 we describe the results of fitting a 
Compton shoulder model and relativistic disk-line 
model respectively, to the excess flux at the base 
of the core of the Fe Ka line. Finally, we present 
our conclusions in §7. 

2. OBSERVATIONS AND DATA 

NGC 3783 was monitored over a period of 
~ 125 days, starting 2001 February 24, with the 
Chandra HETGS during five snapshots. We will 
refer to these snapshots as observations (1) to 
(5). HETGS consists of two grating assemblies, 
a High-Energy Grating (HEG) and a Medium- 
Energy Grating (MEG), and it is the HEG that 
achieves the highest spectral resolution. The 
MEG has only half of the spectral resolution 
of the HEG and less effective area in the Fe-K 
band, so our study will focus on the HEG data. 
The Chandra data were reduced and HEG spec¬ 
tra made exactly as described in Yaqoob et al. 
(2003b). We used only the first orders of the 
grating data (combining the positive and nega¬ 
tive arms). The mean HEG count rates ranged 
from 0.36 — 0.50 ct/s. The exposure times for 
observations (1) to (5) were in the range 165- 
169 ks. The count rate and exposure time for the 
total time-averaged HEG spectrum were 0.36 ct/s 
and 832 ks respectively. Background was not sub¬ 
tracted since it is negligible over the energy range 
of interest (e.g. see Yaqoob et al. 2003a). Note 
that the systematic uncertainty in the HEG wave¬ 
length scale is ~ 433 km s“^ (~ 11 eV) at 6.4 keV 
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Eig. 1 shows the 2-7 keV lightcurves for the five 
observations, showing the count rates summed 
over the —1 and 4-1 orders of both the MEG and 
the HEG, binned at 1024 s. It can be seen that 
the observations, each of duration ~ 2 days, were 
designed to probe timescales of ~ 1 day and less, 
~ 1,2,3, and 4 weeks, 3 months, and 4 months. 
The 2-7 keV continuum variability over the first 
~ 2 weeks of the campaign was confined to the 
flux remaining within ~ ±25% of the mean dur¬ 
ing that period, but in observation (4), about a 
month into the campaign, the overall flux was 
higher than the mean in the first two weeks of 
the campaign by ~ 50% (see Eig. 1). However, 
the variability about the new mean flux during 
observation (4) was still ~ ±30% relative to the 
mean in observation(4). In the final observation 
the 2-7 keV flux began back at its level at the 
beginning of the campaign, but by the end of 
the ~ 2 day observation it had risen by ~ 50%, 
compared to the level at the beginning of final 
observation of the campaign. The excess vari¬ 
ance above the expectation for Poisson noise (e.g. 
Turner et al. 1999), calculated over the entire ~ 4 
month monitoring period from the lightcurve in 
Fig. 1 is (4.52 ±0.28) x lO'^. 


3. SPECTRAL FITTING RESULTS 

We used XSPEC vll.2 (Arnaud 1996) for spec¬ 
tral fitting. Since we were interested in utiliz¬ 
ing the highest possible spectral resolution avail¬ 
able, we used spectra binned at 0.0025%, and this 
amply over-samples the HEG resolution (0.012% 
FWHM). The C-statistic was used for mini¬ 
mization. By definition, calculation of the C- 
statistic requires only knowledge of the number 
of counts in a bin, but for spectral plots, the er¬ 
ror bars shown correspond to asymmetric errors 
calculated using the approximations of Gehrels 
(1986). All model parameters will be referred to 
the source frame, unless otherwise noted. Note 
that since all models were fitted by first folding 
through the instrument response before compar¬ 
ing with the data, the derived model parameters 
do not need to be corrected for instrumental re¬ 
sponse. 


^ http: / / space.mit.edu/CXC / calib/hetgcal.html 
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Figure 1. Lightcurves from the Chandra high-energy grating (HETGS) monitoring campaign of NGC 3783. In this paper, 
the five snapshots are referred to as observations (l)-(5) The first panel shows the first two snapshots, and the remaining 
three panels show one snapshot in each panel. Shown are the 2-7 keV count rates combined from MEG and HEG data, 
utilizing both the —1 and -|-1 orders of the gratings, binned at 1024 s. The zero reference time in this composite lightcurve 
is UT 2001 February 24 19:06:15. The scale on the time axis is such that equal lengths correspond to equal time intervals. 


3.1. Simple Continuum Model 

Our method is first to fit a simple empirical 
model to the continuum and extract parameters 
for the core of the Fe Ka emission line to establish 
whether there is any intra-observation variability 
of the line parameters. Since we are comparing 
line parameters, and the line core is narrow (see 
YP04), this is good enough to guide subsequent, 
more detailed analysis. Accordingly, we will then 
describe more physical (and more complex) mod¬ 
els for the continuum and absorption and show 
the extent of the sensitivity of the absolute pa¬ 
rameters of the Fe Ka emission-line core to details 
of these models. 

Thus, we fitted a simple power law plus Gaus¬ 
sian emission-line model over the 2-7 keV band 
for each of the five spectra. The Gaussian com¬ 
ponent had three free parameters (line center en¬ 
ergy, width, and intensity) so the model had five 
free parameters in total, including the continuum 
slope and normalization. 

Although the Fe Ka line consists of two com¬ 
ponents {Kai and Ka 2 j separated by 13 eV), 
we modeled it as a single Gaussian, since it was 
shown in Yaqoob et al. (2001), that with the 
spectral resolution of the HEG, there is a negli¬ 
gible impact on the measured line width. Some 
broadening may also result from the presence of 
line emission from more than one ionization state 
of Fe. However, we will bear this in mind when in¬ 


terpreting the measured FWHM velocities. Also, 
the use of a single Gaussian (without any attempt 
to model the underlying broad Fe Ka emission) 
has a negligible impact on the measured center 
energy of the core (see §6 and Yaqoob et al. 2001). 
In any case, we will compare the parameters ob¬ 
tained from the simple model with those obtained 
from more complex models (see §3.2, §5, and §6). 

The best-fitting emission-line parameters for 
each spectrum are shown in Table 1 (as well as 
extrapolated 2-10 keV fluxes and luminosities ob¬ 
tained using more complex models for the line 
and continuum, as described in §6). In order that 
the results can used for future statistical analyses, 
statistical errors corresponding to 68% confidence 
for three interesting parameters (AC = 3.506) 
are given. In addition, a more conservative mea¬ 
sure, the 90% confidence, three-parameter range 
(AC = 6.251) for each line parameter is also given 
in Table 1. 

Table 1 shows that the peak energy of the 
Fe Ka emission-line core is tightly constrained 
around 6.40 keV in all five observations, and 
shows no evidence for variability even at the 68% 
confidence level. The line intensity has 68% confi¬ 
dence statistical errors of the order of ~ 30% and 
also shows no sign of variability. The equivalent 
width (EW) however, ranging from ~ 50 — 90 eV, 
does show variability at the 68% confidence level, 
but not at the 90% confidence level. However, 
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any EW variability in this case is simply due to 
the continuum level changing whilst the line in¬ 
tensity remains steady. The 2-10 keV flux and 
luminosity show only ~ 50 — 60% variability over 
the entire monitoring campaign. The most inter¬ 
esting result from Table 1 is that in the first two 
observations the Fe Ka line core appears to be re¬ 
solved at the 90% confidence level, but in the lat¬ 
ter three observations it appears to be unresolved. 
The 90% upper limit on the FWHM in the last 
three observations is not more than 3000 km s“^, 
but for the first two observations it is 4000 km s“^ 
or more. For comparison, the FWHM resolution 
of the HEG at 6.4 keV is ~ 1860 km s“^. 

We report the results of investigating this ap¬ 
parent variability of the width of the Fe Ka line 
core further in §3.3, but first we report the results 
of investigating the dependence of the Fe Ka core 
line parameters on the details of modeling the 
continuum and absorption. 

3.2. Dependence of Fe Ka Line Parameters 
on Continuum and Absorption Model¬ 
ing 

Here we use the ~ 830 ks MEG and HEG 
spectra (i.e. summed over the five observations) 
to investigate in detail the effect of using an 
over-simplified empirical continuum (as described 
above) on the derived Fe Ka core line parame¬ 
ters. It is well known that there is a complex 
photoionized absorber in NGC 3783, with at least 
two components. The absorber has been modeled 
in great detail by several research groups, using 
Chandra and XMM-Newton data (e.g. Kaspi et 
al. 2002; Netzer et al. 2003; Krongold et al. 
2003; Behar et al. 2003). Since highly ionized 
components of the absorber can affect the data 
even in the Fe-K band, it is conceivable that the 
warm absorber can affect interpretation of the 
broad part of the Fe Ka emission line. In fact, 
the apparent lack of a broad Fe Ka emission line 
in recent XMM-Newton data was attributed to 
the complex absorber not modeled adequately in 
previous analyses (Reeves et al. 2004). Reeves et 
al. (2004) also reported an absorption line due 
to He-like or H-like Fe in the same XMM-Newton 
data. We expect the effect of the photoionized ab¬ 
sorber on the narrow Fe Ka emission-line core to 
be unimportant, but since we will also be search¬ 
ing for spectral variability over a broader energy 
range than just the line core, and revisiting the 
question of the presence of a broad Fe Ka line. 
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we constructed a model of the photoionized ab¬ 
sorber. 

We compiled a spectral energy distribution 
(SED) from the literature to use as an input 
to the photoionization code XSTAR^,which was 
then used to construct grids of models for spec¬ 
tral fitting. Full details and results for the com¬ 
plex absorber modeling are given in McKernan, 
Yaqoob, & Reynolds (2005, MN, submitted). For 
the present purpose of Fe Ka line modeling, we 
did not attempt to model the individual soft X- 
ray absorption-lines and features at the level of 
detail accomplished in the dedicated studies men¬ 
tioned above. We found that a two-component 
absorber combined with a broken power-law con¬ 
tinuum was sufficient to model the broadband 
MEG and HEG data. The model was derived 
by first fitting the MEG 0.5-7 keV data, using an 
additional Gaussian to model the Fe Ka emission¬ 
line core, whose parameters were initially frozen 
at the values derived using the simple power- 
law continuum model (§3.1). Then this model 
was fitted to the 0.8-9 keV HEG data with all 
parameters except the hard X-ray slope frozen. 
Then, restricting the energy range to 2-7 keV 
(in order to directly compare derived emission¬ 
line parameters with the empirical continuum fits 
in §3.1) the Gaussian emission-line parameters 
(peak energy, width, and intensity) were allowed 
to float and the best-fit found. The photoion¬ 
ized absorber components had column densities 
of 4.7 X 10^^ cm“^ and 3.5 x 10^^ cm“^, with 
corresponding ionization parameters log^ = 0.74 
and 2.21 respectively (where ^ = Lion/ugr^, Lion 
being the 1-1000 Rydberg ionizing luminosity, Ug 
the electron density, and r the source to absorber 
distance). The soft X-ray modeling is discussed in 
detail in McKernan et al. (2005, MN, submitted), 
as well as in the several studies mentioned above 
so we do not discuss it any further. Fig. 2 (a) 
shows the ratio of the summed HEG data to this 
model in the 2-7 keV band, with the data binned 
at O.OOSA. The other five panels in Fig. 2 (b)-(f) 
show the ratio of the data, from each of the five 
snapshots, to this same model, with only the over¬ 
all normalization adjusted to obtain the best-jit. 
Aside from line-like features between ~ 6 — 7 keV, 
the model gives a strikingly good fit to the spec¬ 
tra from all five of the separate observation data 
sets. We will address the possible origin of these 


^http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.html 
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Figure 2. Ratios of data/model for the Chandra HEG data for NGC 3783. The model consists of a complex continuum plus 
warm absorber and simple Gaussian Fe Kn emission line, fitted as described in §3.2, to the total ~ 830 ks HEG spectrum 
(i.e. all snapshots combined). The ratio of the total HEG spectrum, binned at 0.005A, to this model is shown in (a). The 
ratios of HEG spectra (for the individual snapshots binned at O.OlA) to exactly the same model (with all parameters fixed 
at the values obtained from fitting the total spectrum, apart from an overall renormalization) are shown in (b)-(f). Note 
that the bin size for the individual snapshot spectra in (b)-(f) is comparable to the HEG FWHM spectral resolution of 
~ O.OI 2 A. It can be seen that there is very little evidence for spectral variability, aside from line-like residuals between 
~ 6 — 7 keV in some of the ratio plots. Vertical dashed lines indicate energies of 6.400 keV, 6.700 keV, and 6.966 keV in 
the NGC 3783 rest frame, corresponding to the energies of the Fe I Ka, Fe xxv — ls2p (resonance), and Fe xxvi Ly« 
transitions. The arrow in (a) indicates the position in energy of the Fe-K absorption feature found in the XMM-Nemton 
data by Reeves et al. (2004). 
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residuals in more detail in §6. The apparently 
variable line-like features at ~ 7 keV are likely to 
be due to Fe xxvi Lyct emission (which, blended 
with Fe K/3, was detected in XMM-Newton data 
by Reeves et al. 2004; also see Kaspi et al. 2002). 

Fig. 3 shows joint, two-parameter, 99% con¬ 
fidence contours directly comparing the derived 
Fe Ka line results for the simple and complex 
continuum models. Black contours show the re¬ 
sults using the simple, empirical power-law con¬ 
tinuum, and red contours show the results using 
the warm absorber plus complex intrinsic con¬ 
tinuum described above. Specihcally, Fig. 3 (a) 
shows the Fe Ka line intensity versus line cen¬ 
ter energy contours, whilst Fig. 3 (b) shows the 
EW versus FWHM contours. It can be seen that 


although both the line intensity and EW are con¬ 
sistent at the 99% confidence level from the two 
different continua, the fits using a complex con¬ 
tinuum give values for the best-fitting line inten¬ 
sity and EW which are ~ 20% higher than those 
obtained using the simple power-law continuum 
only. This is because the intrinsic flux from the 
line has to be larger to compensate for the absorp¬ 
tion in the complex model. However, we do not in 
fact know whether the warm absorber lies farther 
from the X-ray continuum source than the Fe Ka 
line emitter so there is an inherent uncertainty in 
the Fe Ka line intrinsic intensity which will only 
be resolved when we are able to constrain the ge¬ 
ometry of the system better with future missions 
(but see §5 and §6). 
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Figure 3. Joint, two-parameter, 99% confidence contours for the Fe Ka emission-line core from the total, 830 ks, Chandra 
HEG spectrum of NGC 3783, illustrating the effects of different levels of approximation in modeling the continuum. Black 
contours correspond to a simple power-law continuum (fitted in the 2-7 keV band), red contours correspond to a model 
including a photoionized absorber, as described in §3.2, and blue contours include a Compton-refiection continuum in 
addition to the photoionized absorber, also described in §3.2. The emission-line is modeled by a simple Gaussian and (a) 
shows the line intensity versus center energy, and (b) shows the line equivalent width (EW) versus FWHM. All quantities 
are in the rest frame of NGC 3783. It can be seen that the center energy of the line is not sensitive to details of modeling 
the continuum, but there is some model dependency on the line intensity, EW, and FWHM. However, the differences are 
not statistically significant at the 99% confidence level. The dotted line in (a) corresponds to the rest-energy of Fe i Aq 
( 6.400 keV) and the dotted line in (b) corresponds to the HEG FWHM resolution at the observed center energy of the 
Fe Ko line. 


The Chandra HEG data are not sensi¬ 
tive enough to strongly constrain a Compton- 
refiection continuum. However, we investigated 
the effect of including a Compton-refiection con¬ 
tinuum (in addition to the complex ionized ab¬ 
sorber model described above) since it is conceiv¬ 
able that it could still affect the derived line pa¬ 
rameters as it has a complex shape in the 6-7 keV 
region due to the Fe-K edge. Moreover, if the 
Fe Ka emission line is formed in optically-thick 
matter (for example reflection from the inner sur¬ 
face of a putative obscuring torus), a Compton- 
refiection continuum is expected. We used the 
PEXRAV model in XSPEC (see Magdziarz and 
Zdziarski 1995). Although this model is based 
on a disk geometry, it is adequate for our pur¬ 


pose since the Chandra data are not sensitive to 
the differences in geometry. Again, since the data 
cannot constrain this model well, the parame¬ 
ters of the reflection model were fixed at nomi¬ 
nal values (except for the photon index and nor¬ 
malization of the intrinsic power-law continuum). 
Specifically, we fixed the Fe abundance at the de¬ 
fault solar value, the disk inclination angle at 30°, 
and the effective solid angle of the reflector at 27r 
(i.e. the “reflection factor”, R = 1). However, in 
§5 and §6 we will investigate the effect of allowing 
R to be free. The e-folding cutoff energy of the 
power law was fixed at 300 keV, well outside the 
range of the data. 

Fig. 3 (a) shows the joint, two-parameter, 99% 
confidence contour (blue) of the Fe Ka line inten- 
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sity versus line center energy, and Fig. 3 (b) shows 
the two-parameter, 99% confidence contour of the 
line EW versus FWHM. These contours are di¬ 
rectly compared with the corresponding contours 
derived with the Compton-reflection continuum 
omitted (red), and with only a simple power-law 
(black). It can be seen that the effect of omitting 
the Compton-reflection continuum on the derived 
parameters of the core of the Fe Ka emission line 
is negligible, even at the 99% confidence level. 
The effect of the ionized absorber on the line pa¬ 
rameters is far more important. 

3.3. Variability of the Fe Ka Line Core 

Fig. 4 (a) shows the joint, two-parameter, 99% 
confidence contours of the Fe Ka line intensity 
versus center energy obtained from the time- 
averaged HEG spectrum (black) compared with 
those obtained from the five individual observa¬ 
tions (or “snapshots”). This confirms the result 
from the spectral fitting (see Table 1), that at 
99% confidence, the intensity and peak energy of 
the Ee Ka line core do not vary. 

The five snapshots span a range of timescales 
from < 2 days to ~ 125 days. However, long-term 
intensive monitoring of NGC 3783 with RXTE, 
some of which was simultaneous with the present 
Chandra observations, shows that the continuum 
varies significantly on all timescales in this range 
(Markowitz & Edelson 2004). Therefore all we 
can say is that the actual response time of the 
Fe Ka line core must be greater than ~ 2 days 
because if the continuum varies more rapidly than 
the response time, the Fe Ka line core intensity 
measured will correspond to some average con¬ 
tinuum level. If we take t > 169 ks as the re¬ 
sponse time for the Fe Ka line core to continuum 
variations, we get r > 5 x 10^® cm for the loca¬ 
tion of the line-emitter relative to the continuum 
source (the recombination timescale likely does 
not play a role here since the center energy of the 
Fe Ka line core indicates neutral Fe). For a cen¬ 
tral black hole mass of 2.9 x 10^ Mq (Peterson et 
al. 2004), this corresponds to r > lOOOrg, where 
Tg = GMjf? is the gravitational radius. However, 
the FWHM of the Fe Ka line core (~ 1700 km s“^ 
from the total combined spectrum) places its ori¬ 
gin much further from the central black hole. The 
actual width already indicates that it might coin¬ 
cide with the outer BLR/inner NLR. Note that in 
the present paper, all quoted line widths derived 
from spectral fitting are already corrected for the 


instrumental line-spread function. Assuming a 
virial relation and an r.m.s. velocity dispersion 
of ■\/3ufwhm/2 (e.g. Netzer 1990), places the line 
emitter at 0.06 pc, or ~ 70 light days from the 
central continuum source (using the above mass). 
If some of the line width is not kinematic in ori¬ 
gin (e.g. if there are multiple lines from differ¬ 
ent ionization states of Fe), then this distance is 
even greater. This is consistent with the lack of 
response of the line intensity to continuum vari¬ 
ability. However, the lack of variability of the nar¬ 
row Fe Ka line core is not necessarily an expected 
result, because variability in the narrow Fe Ka 
line core has been observed in Mkn 841 (Petrucci 
et al. 2002), independent of the continuum (and 
that result is still not fully understood, but see 
Longinotti et al. 2004). 

In §3.1 (Table 1) we showed that observations 

(1) and (2) appeared to have a somewhat broader 
Fe Ka line core than observations (3)-(5). To 
investigate further, we combined the five HEG 
spectra into two groups, one spectrum from ob¬ 
servations (1) and (2), and the other spectrum 
from observations (3)-(5). Fig. 4 (b) shows the 
joint, two-parameter, 68%, 90%, and 99% con¬ 
fidence contours of the Fe Ka line EW versus 
the FWHM, for these two grouped spectra, com¬ 
pared with the 99% contour obtained from the 
mean, time-averaged spectrum. Indeed, the 68% 
confidence contours from observations (1) and (2) 
do not overlap with those from observations (3)- 
(5), the latter showing a narrower line width than 
the former, with the narrower line having a some¬ 
what smaller flux than the broader line. Thus, at 
68% confidence, the Fe Ka line core appears to be 
broader and more intense in observations (1) and 

(2) than in observations (3)-(5). However, the 
90% and 99% confidence contours do overlap so 
we cannot say that the line width varied at 90% 
confidence or greater. On the other hand, the 
line is resolved by the HEG at higher than 99% 
confidence in observations (1) and (2) but is not 
resolved at 99% confidence in observations (3)- 
(5). Future instrumentation will have to address 
the question of variability in the line width. 

4. He-like Fe Absorption Feature 

Kaspi et al. (2002) reported He-like Fe ab¬ 
sorption in the Chandra HETCS data. Reeves 
et al. (2004) detected the absorption line from 
an XMM-Newton observation of NGG 3783 (in 
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Figure 4. (a) Joint, two-parameter, 99% confidence contours of line intensity versus center energy (both in the AGN rest 
frame) for the Fe Ka emission-line core (modeled with a simple Gaussian) from Chandra HEG spectra of NGC 3783. The 
continuum was modeled by a simple power-law, fitted over the 2-7 keV energy range. The small, black contour was obtained 
from the total, 830 ks, HEG spectrum. Remaining contours were obtained from the individual snapshots as follows: black 
(large), observation (1); red, observation (2); green, observation (3); blue, observation (4); magenta, observation (5). See §3 
for details. It can be seen that there is no discernible variability at the 99% confidence level. The dotted line is at 6.400 keV 
in the rest frame of NGC 3783. (b) Joint, two-parameter, confidence contours of line equivalent width (EW) versus line 
FWHM (both in the AGN rest frame) for the Fe Kq emission-line core in the NGC 3783 HEG data, using the same model 
as in (a). Black contour shows the 99% confidence region for the total, 830 ks, HEG spectrum. Remaining contours were 
obtained by combining the snapshots from observations (1) and (2) (blue contour), and observations (3), (4), and (5) (red 
contour). Dotted lines show the 68%, and 90% confidence regions and solid lines show the 99% confidence contours. These 
combinations were formed after examining the contours of all five observations separately and grouping them according to 
similar widths. It can be seen that there is marginal evidence for a change in the effective width of the line, which appears 
to be resolved in observations (3), (4) and (5), but not in observations (1) and (2) (the dotted line corresponds to the HEG 
FWHM resolution at the observed center energy of the line). 

December 2001), centered at 6.67 ± 0.04 keV (in 
the source rest frame), with an EW of 17 ± 5 eV. 

The line was attributed to Fe XXV Is^ — ls2p 
resonance absorption and it was not resolved by 
the CCDs. Evidence for this absorption feature 
can be seen in Fig. 2 (an apparent absorption dip 
is indicated by an arrow). We modeled the fea¬ 
ture by adding an inverted Gaussian to the two- 
component warm absorber plus broken power- 
law model described in §3 (including the Gaus¬ 
sian emission line). All parameters except for the 
hard power-law slope, the overall continuum nor¬ 
malization, and the three absorption-line param¬ 
eters (center energy, intrinsic width) were frozen 


at their best-fitting values. The C-statistic de¬ 
creased by 9.2 upon the addition of the absorp¬ 
tion line, which corresponds to a detection signifi¬ 
cance of only 99% for three additional free param¬ 
eters. For the line center energy, EW, and intrin¬ 
sic width, we measured 6.62lQ;Qg keV, 12^® eV, 
and 6405^2670 ^m s“^ FWHM respectively (er¬ 
rors are one-parameter, 90% confidence, to facil¬ 
itate direct comparison with the XMM-Newton 
results of Reeves et al. 2004). Thus, the Chan¬ 
dra HETGS and XMM-Newton measurements are 
consistent with each other. 
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Figure 5. Ratios of HEG data for NGC 3783 to a simple 
power-law continuum model, fitted in the 5—7 keV band, 
omitting the 6.0-6.5 keV data. Black data points corre¬ 
spond to the summed spectra from observations (1) and 
(2), and the red data points correspond to the summed 
spectra from observations (3), (4), and (5). These two 
groupings show similar line widths within each group, from 
examination of the data from the individual observations 
(see Fig. 4 and §3). This direct comparison shows how 
the Fe Ka line core is slightly broader in observations (1) 
and (2) (black) than in observations (3), (4), and (5) (see 
Fig. 4). The continuum was fitted to the two data sets 
simultaneously, with the photon index tied together. Al¬ 
though the continuum is over-simplified, the purpose here 
is simply to compare the emission-line core for the two 
data sets. 


5. Fe Kof Line Compton Shoulder 

Using the same NGC 3783 Chandra HEG data 
as described in the present paper, Kaspi et al. 
(2002) noted the detection of a so-called “Comp¬ 
ton shoulder” on the red side of the peak of 
the Fe Ko emission line. This structure would 
be due to the core Fe Ka emission-line photons 
Compton-scattering on electrons in the medium 
in which the line is formed, before escaping the 
medium and reaching the observer. If the ma¬ 
terial in which the line is formed is not too hot 
(i.e. kT <C 6.4 keV), the Compton scattering 
results in the line photons losing energy so that 
the scattered emission appears redward of the 
peak of the line. The relative magnitude of the 
Compton shoulder relative to the peak line flux 
increases with optical depth of the medium, un¬ 
til the medium becomes optically-thick (see, for 
example. Matt 2002, for a theoretical description 
and Monte Carlo calculations). The peak of the 
once-scattered photon energy distribution is ex- 
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pected at two Compton wavelengths longer than 
the wavelength of the unscattered line photons. 
For a line at 6.4 keV, this is ~ 6.25 keV in the 
rest frame of NGC 3783, or 6.19 keV in the ob¬ 
served frame. The excess flux to the red side of 
the Fe Ka line core can be seen in Fig. 2 (a) and 
Fig. 5 and does in fact appear to extend down to 
the appropriate energy. 

We applied a model of the Compton shoulder 
to the NGC 3783 HEG data. For the line scatter¬ 
ing we used the XSPEC model of Watanabe et al. 
(2003). The intrinsic width of the emission line 
in this model is much less than the instrument 
resolution and includes no kinematic information 
so is rather unphysical. During preliminary fit¬ 
ting we found that the model always left excess 
flux blueward of the Fe Ka line core. There¬ 
fore we convolved the line profile with a Gaus¬ 
sian whose width was a free parameter, in or¬ 
der to mimic Doppler broadening of the intrinsic 
line profile. We used the broken power-law plus 
two-component warm absorber, with a Compton- 
reflection continuum as the baseline model, as de¬ 
scribed in §3.2. The parameters of the warm ab¬ 
sorber and Compton reflection models were fixed 
at the values described in §3.2. The fitting was 
performed between 2-9 keV (i.e. the high-energy 
end of the fitted range was now extended from 
7 keV up to 9 keV). Only the hard X-ray slope, 
the overall continuum normalization, and the col¬ 
umn density of the Compton-scattering medium 
were free parameters. Since the Fe Ka line peak 
energy indicates low ionization states of Fe and 
is consistent with Fe i, the temperature of the 
scattering medium must be low and we assumed 
the medium is neutral (“cold”). The data, best¬ 
fitting model, and data/model ratio are shown 
in Fig. 6 (a). Note that the He-like Fe ab¬ 
sorption line discussed in §4 has deliberately not 
been modeled in Fig. 6 (a) in order to clearly 
show its presence and magnitude, especially in 
the data/model ratios, and to facilitate compari¬ 
son with the XMM-Newton data. We obtained a 
FWHM of I 5 OOI 34 Q km s“^ for the width of the 
Gaussian used to convolve the Compton-scattered 
line profile. The column density we derived was 
7.5to'6 ^ 10^^ cm“^ (90%, one-parameter errors). 
This corresponds to a Thomson depth of ~ 0.60, 
which is not quite Compton-thick. The Thom¬ 
son depth of the scattering medium must be less 
than ^ or else the mean number of scat¬ 
terings would be ~ 2 and the peak energy of 
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the scattered photons would be much lower than 
measured (the scattered line photons would have 
a distribution with significant flux down to four 
Compton wavelengths below the zeroth-order line 
peak). A Thomson depth of -\/2 corresponds to 
~ f.7 X fO^^ cm“^. 

Ideally, we would self-consistently model the 
Compton-reflection continuum and the Compton- 
scattered line since both are produced by the 
same physical process. When we allowed the rela¬ 
tive normalization of the Compton reflection con¬ 
tinuum, R, to be a free parameter, we obtained 
R = l.OtoJ (justifying i? = 1 in the fits thus 
far). However, the Compton-reflection continuum 
model is calculated for an optically-thick disk, but 
the emission line may be produced in transmis¬ 
sion since we measured a small Thomson depth 
from the Compton shoulder. Better data are re¬ 
quired to justify more sophisticated modeling. 

6. Relativistic Disk Line Model 

We investigated whether a relativistically 
broadened Fe Ka line (in addition to a narrow 
Gaussian component) could also account for the 
Fe Ka line profile in NGC 3783. We fitted the 2- 
9 keV HEG data using a broken power-law plus 
two-component warm absorber, with a Compton- 
reflection continuum as the baseline model, as de¬ 
scribed in §3.2. Initially, the Compton-reflection 
parameter, R, was fixed at unity, which corre¬ 
sponds to the expected steady-state normaliza¬ 
tion of the reflected continuum from a neutral 
(“cold”) Compton-thick disk subtending a solid 
angle of 27r at the X-ray source. Preliminary 
spectral fitting with R free showed that the best¬ 
fitting value is i? ~ 1 in any case (as in §5), and 
we will give statistical errors when R was a free 
parameter below. In addition to the above model 
components we included a Gaussian (to model the 
line core), and an emission line from a relativistic 
disk around a Schwarzschild black hole (e.g. see 
Fabian et al. 1989). The energy, intrinsic width, 
and intensity of the Gaussian were free parame¬ 
ters. The inner radius of line emission from the 
disk was initially fixed at (r^ = GM/c^) and 
the outer radius was fixed at lOOOrg. The line 
radial emissivity was assumed to be a power law 
(line emissivity proportional to with the in¬ 
dex, < 7 , free. In the disk rest-frame the emission¬ 
line energy was fixed at 6.4 keV. The disk incli¬ 
nation angle and overall line intensity were also 
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free parameters. Thus, there were a total of seven 
free parameters for this model, including the con¬ 
tinuum slope and normalization. 

The data, best-fitting model, and data/model 
ratio are shown in Fig. 6(b). It can be seen 
that the fit is as good as that for the Compton¬ 
scattering model, shown in Fig. 6(a). Again, we 
deliberately did not model the He-like Fe absorp¬ 
tion feature discussed in §4 in order to clearly 
show its effect and to facilitate comparisons be¬ 
tween Chandra and XMM-Newton data. In the 
5.7-6.9 keV energy band shown in the figure, the 
C-statistic for the disk-line model is less than 
that for the Compton shoulder model (§5) by 
only 1.6 for an additional three free parameters. 
Thus, although the disk-line model is not sta¬ 
tistically preferred over the Compton shoulder 
model, it cannot be ruled out. The center en¬ 
ergy of the Gaussian was 6.399 keV, with statis¬ 
tical errors similar to those obtained from sim¬ 
pler continuum models (Table 1). Since the peak 
energy was well-determined, this parameter was 
fixed when deriving the statistical errors for the 
other parameters in the model (in this section, 
all errors are 90% confidence, one-parameter, or 
AC = 2.706). The Gaussian line FWHM ob¬ 
tained was 1460lggQ km s“^. The Gaussian line 
width was then fixed at the best-fitting value in 
order to derive statistical errors on the remaining 
model parameters, in order to avoid the spectral 
fitting becoming unstable. 

We obtained EWs of 56tii eV and eV for 

the Gaussian and disk emission-line components 
respectively (all line measurements are in the 
source frame). The radial emissivity index mea¬ 
sured for the disk-line emission was q = I.fiOto^gr, 
and the inclination angle was constrained with an 
upper limit only, of 11°. When we allowed the 
relative normalization of the Compton-reflection 
continuum, R, to be a free parameter we ob¬ 
tained R = l.Oto^a- recall that i? = 1 

corresponds to the case of time-steady illumina¬ 
tion of a thick disk subtending 27r solid angle 
at the X-ray source. We note that De Rosa et 
al. (2002) obtained R = 0.71i'ig’28 from non- 
contemporaneous BeppoSAX data for NGC 3783, 
and Markowitz, Edelson, & Vaughan (2003) ob¬ 
tained R = 0.62 ± 0.12 from RXTE data which 
overlapped with the first three Chandra observa¬ 
tions reported in the present paper, consistent, 
within the errors, with the Chandra value. How¬ 
ever, Markowitz et al. (2003) also reported R = 
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Figure 6. (a) Modeling of the Compton shoulder of the Fe Ka line in the total, 830 ks, Chandra HEG spectrum 

of NGC 3783 (the data are binned at 0.005 A for the plot). The model consists of a power law, photoionized absorber, a 
Compton-reflection continuum, and a narrow Fe Kq line which suffers Compton scattering in cold, solar-abundance material, 
observed in transmission. The line profile was convolved with a Gaussian, whose best-fitting width was ISOO^^^q km s“^. 
See §5 for details of the fitting procedure and results, (b) The Chandra HEG spectrum, same as in (a), fitted with an 
emission line profile from a relativistic accretion disk rotating about a Schwarzschild black hole, plus a Gaussian to model 
the Fe Ka line core. See §6 for model-fitting details and results. Note also that the dip at 6.51 keV (observed frame) is 
not real: it is narrower than the instrument resolution (0.012A FWHM) and is not detected in boththe —1 and -hi orders of 
the grating, (c) The Chandra Compton-shoulder model in (a) fitted to non-contemporaneous XMM-Newton EPIG pn GCD 
data for NGC 3783, with only the overall normalization and power-law slope free, (d) The Chandra disk line plus Gaussian 
model in (b) fitted to non-contemporaneous XMM-Newton pn data for NGC 3783, with only the overall normalization and 
power-law slope free. See Reeves et al. (2004) for details of the XMM-Newton data. We note that all of the line profiles 
shown here are not unfolded: the profiles correspond to the ratio of counts to predicted model counts, multiplied by the 
best-fitting model. Note that the absorption line at ~ 6.7 keV was deliberately not modeled for any of the four cases shown 
here, in order that the reader can readily compare the absorption line in the Chandra and XMM-Newton data. See §4 or 
Reeves et al. (2004) for detailed fitting of the absorption line in the Chandra or XMM-Newton data respectively. 


0 .47^0 04 RXTE data taken over a ~ 3.2 year 
period between 1999-2002. Over such a long 
timescale variability can obviously be a factor re¬ 
sponsible for the smaller value of R. 

Upon allowing the outer disk radius to be a 
free parameter in the original disk model, with 
R fixed at 1.0, we obtained a lower limit on its 
value of ^ 540rg. The inner disk radius and q are 
correlated in the sense that they affect the overall 
width of the emission line: larger inner radii allow 
steeper radial emissivity index (g) since the part 
of the disk giving the broadest part of the line 
with a steep q would then be missing. The data 
do not allow the ambiguity between q and the 


inner radius to be removed. By examining con¬ 
fidence contours of q versus rj^, we found that a 
radial emissivity steeper than even r~^ is allowed 
if Hn is of the order of lOOrg or larger. 

We note that our results for the relativistic 
disk line are not inconsistent with the findings 
of Reeves et al. (2004), namely that the XMM- 
Newton data do not require a disk line when the 
warm absorber is taken into account. Reeves et 
al. (2004) obtained a 90% upper limit on the EW 
of a disk line of 35 eV. However, this assumed a 
steep emissivity index (g = 3) combined with a 
small inner disk radius (rin = 6rg), and a disk 
inclination angle of 30°. With a flatter emissivity 
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and a smaller inclination angle, as obtained from 
the Chandra HEG data, a larger EW is allowed 
because the disk line width is narrower. Thus, 
the Chandra HEG and XMM-Newton measure¬ 
ments are consistent with each other. The reason 
why the XMM-Newton data do not yield a signif¬ 
icant detection of the disk line is that the Chan¬ 
dra HEG data are more sensitive to a narrow disk 
line, due to the factor of ~ 4 better spectral reso¬ 
lution. In Fig. 6 (c) and Fig. 6 (d), we show that 
both the best-fitting Chandra HEG Gompton- 
shoulder model and disk-line model respectively, 
fit the XMM-Newton data equally well. For fit¬ 
ting the XMM-Newton data (the same data as in 
Reeves et al. 2004, which were taken at a differ¬ 
ent time to the Chandra data), all parameters ex¬ 
cept the overall normalization and the power-law 
slope were fixed at their Chandra best-fitting val¬ 
ues. Again, the He-like Fe absorption feature dis¬ 
cussed in §4 was deliberately not modeled here in 
order to facilitate comparison between the Chan¬ 
dra and XMM-Newton data. 

Although Fe K(3 line emission is expected to 
accompany Fe Ka line emission, the HEG data 
are not so sensitive to Fe K(3 lines because the 
branching ratio is 17:150 for K(3 : Ka. Having 
said that, line-like residuals are apparent in the 
spectral ratios shown in Fig. 2 but some part of 
this line emission could be due to Fe xxvi Lya. 
Reeves et al. (2004) reported significant line 
emission at ~ 7 keV from XMM-Newton data and 
obtained an EW of 20 ± 5 eV for Fe xxvi Lya 
after accounting for Fe K/3 emission. Since Fe K(3 
line emission is expected to accompany the Fe Ka 
line emission, we self-consistently added Gaussian 
and disk-line Fe K/3 emission with the expected 
branching ratio to the Gaussian plus disk-line 
model described above in order to test the Chan¬ 
dra HEG data for Fe xxvi Lya emission. We 
modeled the latter with an additional Gaussian, 
fixing the energy at 6.966 keV (e.g. see Pike et al. 
1996) but leaving the intrinsic width free for the 
summed, ~ 830 ks HEG spectrum, but fixing it at 
the resulting best-fit value for the individual spec¬ 
tra from observations (l)-(5). The best-fitting 
intrinsic width from the summed HEG spectrum 
was 9200 km s“^ FWHM. However, we did not 
obtain a very significant detection (AC was al¬ 
ways less than 2.7) for either the summed spec¬ 
trum (90%, one-parameter upper limit: 17 eV) or 
the individual spectra from the five observations. 
Although Fig. 2 shows evidence of variability of 
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the line-like residuals at ~ 7 keV, we could only 
obtain upper limits on the EW (29, 27,65,41, and 
14 eV for observations (1), (2), (3), (4) and (5) 
respectively). All of the values of the Fe xxvi 
Lya EW upper limits we obtained are consistent 
with the EW measurements from XMM-Newton 
data by Reeves et al. (2004). 

7. CONCLUSIONS 

We have presented results of X-ray spec¬ 
troscopy of the Fe-K line region in NGC 3783 
from an extended monitoring campaign with 
the Chandra High Energy Grating Transmission 
Spectrometer {HETGS). Consistent with previ¬ 
ous studies, the Fe Ka line core is resolved in 
the time-averaged spectrum by the High Energy 
Grating (HEG) and has a center energy indi¬ 
cating an origin in neutral or lowly-ionized Fe. 
Despite a factor of ~ 1.5 variation in the X- 
ray continuum luminosity, we measured no vari¬ 
ability in the intensity of the Fe Ka line core 
during five observations comprising the observ¬ 
ing campaign, over ~ 125 days. The lack of re¬ 
sponse to the continuum is consistent with the 
fact that a virial interpretation of the line FWHM 
(~ 1700 km s“^) gives a distance between the pu¬ 
tative central black hole and the line emitter of 
at least 0.06pc, or 70 light days. However, we 
detected marginal evidence for variability in the 
intrinsic width of the Fe Ka line core. In two out 
of the five snapshots the Fe Ka line was resolved 
and had a FWHM of ~ 2300 km s“^, whilst in 
the other three snapshots the line was unresolved. 
This apparent variability was not related to the 
X-ray continuum level. At 68% confidence, the 
Fe Ka line core was broader and had a higher 
flux in the first two ~ 170 ks snapshots than in 
the last three. If the line width variability is real 
it could be due to changes in the matter distri¬ 
bution and/or ionization state as a function of 
distance from the central engine and would be 
important to investigate with future missions. 

We found that the excess flux around the base 
of the Fe Ka line core can be modeled by an 
emission line from a relativistic disk rotating 
around a black hole, as well as a model in which 
line core photons are Gompton-scattered (form¬ 
ing a “Compton shoulder” on the red side of the 
line peak). More realistic modeling should self- 
consistently include Compton-scattering of line 
photons, whether the line photons originate in an 
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optically-thick disk or more distant matter. Also, 
the Compton-scattered continuum should be self- 
consistently computed. In principle, measure¬ 
ment of the unscattered and scattered line pho¬ 
tons, along with the reflected continuum, would 
constrain whether the core of the Fe Ka line orig¬ 
inates in optically-thick matter (such as the puta¬ 
tive obscuring torus) or in optically-thin matter 
(e.g. the BLR and/or NLR). However, the signal- 
to-noise of the current data is not sufficient to dis¬ 
tinguish between these scenarios. We note that a 
Compton shoulder has not been detected in any 
other Chandra HETGS observation of a Seyfert 1 
galaxy (Yaqoob & Padmanabahn 2004). If this 
is verified by higher signal-to-noise data, it would 
suggest that the Fe Ka line core is in general likely 
to originate in optically-thin matter. In the case 
of NGC 3783, the EW of the Fe Ka line core de¬ 
duced from the composite Gaussian plus disk-line 
model can easily be produced by optically-thin 
matter with Nh ~ 10^^ cm“^, covering less than 
half the sky (e.g. see Yaqoob et al. 2001). 

Modeled with a relativistic disk line, the disk 
inclination angle was constrained to be less than 
11 ° and the emission line flux from the disk is 
comparable to that in the Fe Ka line core. The 
disk-line modeling included a continuum with a 
complex photoionized absorber and Compton re¬ 
flection modifying the intrinsic continuum. We 
showed that the Chandra HEG data are consis¬ 
tent with the XMM-Newton data in the sense 
that the same disk-line model can account for 
both data sets, even with these continuum com¬ 
plexities. Of course, with the current data it is 
not possible to rule out more complicated sce¬ 
narios involving, for example, a disk line modi¬ 
fied by Gompton scattering, as mentioned above. 
We also confirmed previous detections {HETGS: 
Kaspi et al. 2002; XMM-Newton: Reeves et al. 
2004) of an absorption feature in the Fe-K band 
(at 6.67±0.04 keV, and with an EW of 17±5 eV), 
likely to be due to He-like Fe absorption. On the 
other hand, evidence for Fe xxvi Lya emission 
was marginal in the HEG data. 

Future observations with higher spectral reso¬ 
lution, as afforded by Astro-E2, will be able to 
resolve important ambiguities remaining in the 
Fe-K band in NGC 3783. In particular, it will 
be possible to determine if the Fe Ka line core is 
composed of more than one line component and 
to determine if the broad part of the Fe Ka line 
is really due to a relativistic disk. 
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